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We derive the quantum dynamics of heavy quark systems in quark-gluon plasma as open 
quantum systems. The scatterings of a heavy quark with the hard medium particles give 
rise to Debye screened QCD force and drag force and its fluctuation. We present a 
(^J , unified quantum description of these in-medium QCD forces at high temperature in the 

^ ■ leading-order perturbation on the basis of the influence functional formalism. 
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1. Introduction 

r*** | Heavy quarks play important roles in probing properties of the quark-gluon plasma 



(QGP). A single heavy quark in the QGP interacts with the medium particles 
and eventually it will be kinetically thermalized as is described by the Langevin 
equation!^ 



CO 

o _ 

ft = ~2MT^ + ?®> = ~ (1) 

If one can observe the degree of thermalization in a given time scale, one can derive 
the drag force of the heavy quark. When put in the QGP, a quarkonium, which is 
a bound state of a heavy quark-antiquark pair, gets unstable at sufficiently high 
temperature because the plasma constituents screen the color charges and turn the 
confining potential into the Debye screened potentialP 

i^M=(-Yl + V(r;T)^(r,t), V(r;T) = -C F £-e^ . (2) 

The Debye screening in the QGP and the related issues are still under active theo- 
retical and conceptual developments If one can observe the "melting" of quarko- 
nium, one can infer the temperature provided that there is a reliable theory for the 
melting. 

These heavy quark probes are fundamentally different from each other: The drag 
force is a classical description for the non-potential irreversible force while the melt- 
ing is a quantum phenomenon of bound-states in the screened potential. Moreover 
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it is natural to expect that each heavy quark in a quarkonium experiences the irre- 
versible drag process. In a theoretical point of view, it is desirable to have a unified 
quantum description for both of them, whose classical limit yields the Langevin 
dynamics in the single heavy quark case. In other words, we pursue the unified 
quantum description for in-medium forces in the QGP. The theoretical framework 
needed for this purpose is that of open quantum systems. 

The basics of the open quantum systems^ can be summarized as follows. When 
Hilbert space consists of direct product of two Hilbert spaces Utot — %sys ( 8 ) %cnv J the 
Hamiltonian for the total system is in general H to t — H sys <giI cnv +I S y S (3H onv +Hi n t 
and the density matrix of the total system evolves according to the von-Neumann 
equation: 

j t Ptot(.t) = -i{H tot , Ptot (t)}. (3) 

The dynamics of the open quantum system is described by the master equation for 
the reduced density matrix: 

foPredit) = £Pred(*)> Prcd (*) = Tr c„v [p t ot (*)] • ( 4 ) 

In the case of the heavy quarks in the QGP, the system degrees of freedom are the 
heavy quarks and the environmental degrees of freedom are the light quarks and the 
gluons. In this article, we review the perturbative derivation of the master equation 
for the heavy quark systems at high temperature.^ 



2. Influence Functional of Heavy Quarks 

The description of heavy quark systems as open quantum systems can be derived 
on the basis of the closed-time path formalism of non-equilibrium quantum field 
theory.^ Applied to QCD, the partition function is 

Z[rj,fj] = / V[^,iJ},q,q,A]{^,q,A (s = -oo)\p tot \ip,q,A (s = oo)} 



x exp 



c 



ds / d 3 x{C Q cD{ip,ip, q, q,A) + rjip + iprj} 



(5) 



where ip and q are the heavy and light quark Dirac fields, A is the gluon field, and 
gauge fixing is implied in the QCD Lagrangian £qcd- Here C denotes the closed-time 
path contour which starts from t = — oo, extends to t — oo, and then returns to the 
original t = — oo. The initial density matrix (in the Schrodinger picture) is defined as 
p tot = p*^ v ®p sys : At t = — 00, a system in arbitrary state p sys starts to interact with 
the environmental degrees of freedom in thermal equilibrium. States \tp,q,A) and 
(ip, g, A\ are the coherent states, satisfying for example xp(x)\tp, q, A) = ip(x)\?p, q, A). 

The dynamics of the open quantum system is described by the influence func- 
tional obtained after integrating out the environmental degrees of freedom. The 
influence functional SfvIJi^ , j^} i s expressed as a functional of heavy quark color 
currents on the forward time axis C± and j 2 on the backward time axis Ci 
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where C = C\ + C2 , and is expanded in terms of the currents with coefficients being 
the correlation functions of the gluons in thermal equilibrium. In the leading-order 
perturbation and in the non-relativistic limit, S'fv (or iSpy 1 ^ 11 ) consists only of four- 
Fermi terms whose coefficients are the leading-order gluon two-point function G F b 00 . 
The long-time dynamics of heavy quarks is summarized in the following effective 
action Si+2: 



:NR 



c-LONR 
Dp 



= '•'kin "I" ^FV 

= ( d*x [cZi^i) - 



(0) 
(7) 



cLONR 



dtd s xd s y{jf{t,x),j^{t,x)) 



V{x-y) -iD(x-y) 



jf{t,y) 



/ dtd*xd*y W * D %-® 



iD(x~y) -V*(x-y) 
(/f,NR(*, tyffr V) + jf(t, *)j2,NR.(0 y)) , (8) 



where S ab V(x - y) = -ig 2 G F abfl0 (ui = 0,x - y), D(x - y*) = lmV(x - y), and 
denotes the non-relativistic limit of the heavy quark color current. With this 
effective action S'l+2, the partition function is given by 



Z[r)i,r)x,ri2,ri2\= J £>[0>i,0>i, ip 2 ,ip2](ij}i(-oo)\p s \ip2(-oo)) (9) 
x cxp iS 1+2 + i / d 4 x (771-01 + 0i77i - r7 2 02 - 02?72) • 



From Eq. (j9|), the time evolution of (ip\\p s (t)\ip2) is obtained as a "func- 
tional Schrodinger equation" with an effective Hamiltonian ifi+2 derived from 
S'l+2- However, the physically relevant matrix element of the density matrix is not 
WlPsWIV^) but ps t,Q* 1{c) ,Q 2{c) = (Q* 1{c) \p s (t)\Q 2( C )), where Q (c) s are the non- 
relativistic Pauli spinors for the heavy (anti)quarks. It needs to take following techni- 



cal procedures to obtain the functional Schrodinger equation for ps 



^Qi( c ) I Q2( c ) 



In these procedures, introduction of conjugate fields ip2 
metric notation between -01 and -02 • 



0^2 on C2 allows a sym- 



(i) The operator ordering in if 1+2 is determined according to the time ordering of 
the fields in the path integral formula. 



(ii) The functional Schrodinger equation for ps 
Hi + 2 with substitution 



t, Ql( c ); Q2(c) 



is obtained from 



Qh 



SQl 



(<=). 



(Q2(c) 



Q2(c)) 



Q 



2(c). 



.(10) 



'2(c), 



(iii) The above substitution gives ultraviolet divergent contributions in the four- 
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Fermi terms arising from the Coulomb potential at the origin, e.g. 



/ 



d 3 xd 3 yV{x - y) 



I 



d 3 xC F V(6) 



Q*i{c) 



6 



Q\ {c) m a 



(ii) 



SQ* 1{c) (x)_ 

which needs to be renormalized. 
(iv) Instead of (ii) and (iii) , it is equivalent to derive first the renormalized effective 
Hamiltonian H\™ 2 , in which the operators are reordered in the normal order 
and the ultraviolet divergences are renormalized, and then do the substitution 
(ii) to obtain a functional differential operator H[ e + 2 Qi^yQ^M ■ 

We adopt (iv) and obtain the following renormalized effective Hamiltonian -ff i™ 2 : 



rrrcn 
-"1+2 



J d 3 x oM (0 + 1(c) Qi(c)) + Q\ (c) (-^) Qkc 



(12) 



+ 



d 3 x 



a*M 



(Q2(c)Q 



2(c) 



+ Q2(c) ( -— 1 Q 



2M 



'2(c) 



+ - / d 3 xd 3 yN{(jf(x),jf(x)) 



V(x-y) -iD(x-y) 
-iD{x-y) -V*(x-y) 



I 



+ / d 3 xd 3 y 



V x D(x-y) 
AT 



N {jYnr^V) +3f(x)T 2 , m (y)} 



jf(y) 
j?(y) 



where a = 1 + C F V r me d/2M, V mcd = V(0;T) - V(0;T = 0), and N denotes taking 
the normal-ordered product. From J?i™ 2 i tne functional Schrodinger equation is 
derived as 



d 



t:Ql(c)'Q2(c) — Hl+2[Q*(c)i Q 2(c)] PS t,Ql(c)iQ2(c) 



(13) 



The calculation of the gluon two-point functions Gq ob (w, k) at leading order in the 
hard thermal loop resummed perturbation theory yields 



D(r) = Im [V{r)\ = 
Cv g 2 



- 2iTu 



dk 



sin kr 



r(k 2 



,2 \2 



2tt 



dk 



sin k'i 



r(k 2 +cal) 2 ' 



a = 1 



(-WD-tT), ujl=g 2 T 2 



2M 4tt 



3 6 



(14) 
(15) 
(16) 



At this order of perturbation, the effect of t-channel scattering between a heavy 
quark and hard medium particles is taken into account. This interaction contributes 
to the diffusion process of the heavy quarks as well as the Debye screening of the 
color charges. The former is given by the imaginary part of V{r) and D{r) while 
the latter is by the real part of V(r), and thus these in- medium forces are both 
contained in the renormalized effective Hamiltonian H\ c + 2 - 
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3. Real-time Dynamics of Heavy Quarks 



The master equation for the heavy quark systems is obtained from the functional 
Schrodinger equation Eq. (I13[) . The matrix element of ps(t) between the heavy 
quark coherent states (<2*( c )| and \Q 2 ^} * s gi ven as a functional of the heavy quark 
fields by 



PS 



t, Ql(c)'Q2(c) 



<Q; ( c)Ipb(*)IQ3( C ))> 



\ exp 



- d 



|Q* (c) ) = exp 



- d 



•{0(c)(2)o; ( c)(s)} 



(17) 
(18) 

(19) 



Here |n) is the heavy quark vacuum state. By functionally differentiating 



t, Qu c yQ 



(c)>^2(c) 



we obtain, for example, the reduced density matrix for a sin- 



gle heavy quark system: 

p%{t,x,y) = (x,i\p Q (t)\y,j) cx (il\Q\x)p s (t)Q^ (y)\il) 

= ~5Q\*(x)5Qi*($) PS [*' QI W'^W 
and that for a heavy quark- ant iquark system: 
PQQc i ^2 ! yi j jte ) = {xi,i;x 2 , j\p Q {t)\y 1 ,k;y 2 ,l) 

cx (n | Q l (f i ) Q^ c (af a ) p s (t ) (y 2 ) Q fc t (y , ) | fi) 
c5 5 (5 <5 

PS 



,=0 



(20) 



(21) 



*Qi*(*i) 6Q{* e (x 2 ) 6Q l * c {y 2 ) 8Qfc%) 



t, Ql(c)>Qs(c) 



•1(c) -»2(c) 



=0 



The master equation for each of these reduced density matrices is therefore obtained 



from the functional Schrodinger equation for ps t, Q*r c y Q 2 ( c ) m ^q. vM - 

For example, let us examine the real-time dynamics of a single heavy quark 
system at high temperature. The master equation for a single heavy quark system 
reads 



d 



ifaPqfoxiV) = <)(«- a *) M 



V? 



2M 



(22) 



1 f, , S ik Sji 



—id(x — y) + 



Vxdjx - y) Vi - Vy 



P Q l (t,x,y), 



4T 



iM 



and that for the color-traced reduced density matrix pq (t, x, y) = pg (t, x, y) is given 
by 



d 

i-QlPQ{t,x,y) = < (a- a)M + 



V 2 - V 2 

x y 



PQ(t,x,y) 



+C F < ~id(x -y) + 



2M 

V x d(x~y) V x ~V y 



4T 



iM 



p Q (t,x,y). (23) 
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Interestingly, taking short distance limit d(x) s» d(0) + S 2 \0)x 2 /2 in the master 
equation Eq. ([23| reduces it to the well-known master equation of the Caldeira- 
Leggett model for the quantum Brownian motion.^ Using the master equation 
Eq. ([23)) . the following Ehrenfest relations for averaged expectation values (0)(t) = 
Tr, { PQ {t)0} = Jd 3 xd 3 y{x\p Q {t)\y){y\0\x) are derived: 

= ^) = __?_<£), (24) 
dV 1 M rfr 7 2MT^ h y ' 



To be explicit, 



dt\2M I MT\\2M/ 2 



^d(x)\ x=Q = _^l!v 2 G> , 00 (f)U =0 
C F g{T) 2 f d 3 k 



(25) 



9 / (2^3^oo(- = 0,fc), (26) 

which is consistent with the leading-order perturbative calculation of the drag 
forced In this way, we can derive the consequences of the classical Langevin dy- 
namics through the quantum Ehrenfest relations. 



4. Summary 

In this article, we reviewed the derivation of the master equation for the reduced 
density matrix of heavy quark systems at high temperature. The derivation is based 
on the influence functional formalism applied to the heavy quark systems as the 
open quantum systems and is explicitly shown in the leading-order perturbation. 
The essential information is contained in the renormalized effective Hamiltonian 
HY+2 j horn which the master equation is derived. 

Note that this method is applicable to systems with any finite number of heavy 
quarks and is not limited to the leading-order perturbation; once the renormalized 
effective Hamiltonian is derived, it is straight forward to obtain the master equation. 
As explained in Ref.Pthe forward propagator of heavy quarks can also be derived 
from the renormalized effective Hamiltonian. The complex potential between heavy 
quark-antiquark pair is defined in terms of the forward propagator and thus can 
also be derived from our approach to yield consistent result with RefsP^ in the 
leading-order perturbation. 

Here we list several next steps of this approach: 

• Derivation of the renormalized effective Hamiltonian in the next-to-leading or- 
der of perturbation. We expect the next-to-leading order contribution will in- 
troduce new processes, namely the so-called gluo-dissociation processes! 14 1 15 1 

• Consideration of how to extend this formalism to the general non-perturbative 
situation. One of the challenges is to understand how to describe heavy-light 
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meson bound states solely in terms of heavy quark fields, or to examine what is 
the correct choice of the system variables as the open quantum system in such 
cases. 

• Phenomenological application to the dimuon spectra from bottomonia decays 
at the LHCPlf one has a reliable theory of melting, one can study the physics 
of deconfinement by comparing with the experimental data. 
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